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Abstract 



To study the manifestation of the Aharonov-Bohm effect in many-body 
systems we consider the statistical mechanics of the Gross-Neveu model on a 
ring (1+1 dimensions) and on a cylinder (2+1 dimensions) with a thin solenoid 
coinciding with the axis. For such systems with a non-trivial magnetic flux 
(6) many thermodynamical observables, such as the order parameter, induced 
current and virial coefficients, display periodic but non-analytic dependence 
on 8. In the 2+1 dimensional case we further find that there is an interval of 
9 £ (1/3, 2/3) (modulo integers) where parity is always spontaneously broken, 
independent of the circumference of the cylinder. We show that the mean-field 
character of the phase transitions is preserved to the leading order in 1/iV, 
by verifying the ^-independence of all the critical exponents. The precise 
nature of the quasi-particle, locally fermion-like and globally anyon-like, is 
illuminated through the calculation of the equal-time commutator and the 
decomposition of the propagator into a sum over paths classified by winding 
numbers. 
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I. INTRODUCTION 



Due to their direct relevance to the quantum Hall effect |T|, possible connection with high 
temperature superconductivity as well as their intellectual challenges, anyon systems have 
attracted enormous attention in the past several years. Almost exclusively, the formulation 
of anyons is based on the idea of the Chern-Simons construction. In this formulation anyons 
are locally so-called charge-flux tube composites. Physically, the essence involved is nothing 
but the famous Aharonov-Bohm effect ||. It is the interference of the gauge field, (which 
would be a pure gauge when totally isolated), associated with each particle that modifies 
the statistics of the particle such that it interpolates between the boson and fermion lim- 
its. While it is conceptually more satisfactory, the involvement of the Chern-Simons term, 
because of its technical complexity, makes the understanding of anyon properties rather 
difficult, especially in many-anyon systems. 

In this paper we study a system which in many ways captures the main features of 
anyons, while avoiding the inclusion of the Chern-Simons term. In a topologically simple 
configuration space it would be necessary to include the Chern-Simons term to implement 
the fractional statistics. However, when the configuration space becomes topologically non- 
trivial, there is no need to add terms other than the minimal coupling between the gauge 
and the matter fields in order to achieve the desired interference. Thus with the slight 
complication of introducing a topologically non-trivial configuration space, the procedure 
of setting up a system which displays anyon-like behavior is greatly simplified and one can 
study typical effects such as the periodic but non-analytic dependence of thermodynamical 
observables on the statistics parameter. It is this technical simplification which allowed 
us to carry out our study of the thermodynamics of such a system rather explicitly, using 
conventional methods. We hope that our results will shed some light on the thermodynamics 
of Chern-Simons anyon systems. 

The specific model we consider is the Gross-Neveu model M\ with an external constant 
gauge potential defined by 

C = i>[id, - a M (0)]yV + |^) 2 , (1) 

on a ring (1+1 dimensions) and on a cylinder (2+1 dimensions) with ^-direction compactified 
(of length L) to the leading order in 1/JV. The ip-fie\d is a two-component Dirac spinor 
and implicitly has N components in the internal space. The external gauge field a^{9) is 
generated by a thin solenoid of magnetic flux 2tt9 which coincides with the axis of the ring 
or cylinder, that is, a\{6) = 2-kO/L and all other components vanish. Since the Lagrangian 
is fermionic when 8 = 0, the anti-periodic boundary condition in x-direction, if>(x + L) = 
—if>(x), is naturally imposed. Without losing generality we only need to consider 6 £ [0, 1), 
for the integral part of 9 can be safely gauge-transformed away. The periodic dependence 
of physical quantities on 6 is obvious. 

The geometric setting here is such that the particles move in a multi-connected config- 
uration space. Due to the presence of the magnetic flux, we anticipate that the Aharonov- 
Bohm effect will have a profound impact on the statistical mechanics of the system. 
It is possible to formally remove the gauge field in Eq.(|l]) via a gauge transformation 
if) — > exp(i27i9x / L)ip = if)'. However, the expense of eliminating the gauge field is that 
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the new field ip' becomes multi- valued (when 9 ^ or 1/2). The magnetic flux can be 
regarded as an external means of varying the boundary condition. It should be empha- 
sized that, like in the original Aharonov-Bohm experiment, the magnetic field on the ring 
or cylinder is zero. Therefore, the effect of 9 ^ is purely quantum mechanical in nature, 
taking the advantage of L being finite. In the limit L — > oo all the ^-dependence of physical 
observables should disappear, since then the effect of the boundary condition becomes irrel- 
evant. It is also easy to see that single particle states have smooth ^-dependence and that 
the non-analyticity of physical observables arises only from many-body effects. 

Before we dwell into details, let us comment on the choice of the specific microscopic 
Lagrangian, Eq.(|I]). In our opinion, the precise underlying dynamics is not too crucial, as 
long as the qualitative features, such as the existence of a phase transition, are present. 
Therefore, the choice is simply dictated by convenience. However, it should be pointed 
out that the Gross- Neveu type of Lagrangians do have some relevance to certain condensed 
matter systems or models. For example, the Gross-Neveu model in 1+1 dimensions is closely 
related to the continuum limit of the one dimensional Hubbard model at half-filling and the 
quantum spin-one-half antiferromagnetic Heisenberg chain ||, while in 2+1 dimensions it 
models the continuum limit of the so-called chiral spin liquid |j) . Although we do not directly 
address questions on how to make observations in some specifically designed experiments, 
one should always keep these issues in mind. 

The Lagrangian in Eq.fll]) possesses discrete symmetries. In the 1+1 dimensional case, 
the symmetry is the chiral symmetry 

i/>(x,t) e^Vfo *) , (2) 

whereas in the 2+1 dimensional case, the relevant symmetry is parity 

y,t) -»• i^ 2 ip(x, -y, t) . (3) 

Since the operator ipip changes sign under either Eq.(Q) or Eq.(|3|), the order parameters are 
the ensemble averages {fiif)). In other words, when a mass term is generated, the discrete 
symmetry, the chiral symmetry in 1+1 dimensions or parity in 2+1 dimensions, is dynam- 
ically broken. Because of the discrete nature of the symmetry involved, there will be no 
Goldstone mode associated with the dynamical symmetry breaking in our model. For con- 
venience, we introduce a scalar auxiliary field a, which interpolates the composite operator 
—g 2 ijjif}. Then the Lagrangian in Eq.(|I|) is equivalent to 

C = #0„ - a^jyV - - <rW> ■ (4) 

Note that the a-field obeys periodic boundary conditions for any value of 9, though it is odd 
under the transformations Eq.(^) or Eq.(|3]). 

The rest of the paper is organized as follows. In the next section we present the result in 
1+1 dimensions, along with its relevance to Liischer's small- volume expansion in asymptot- 
ically free field theories with phase transitions. In section III we present the result in 2+1 
dimensions, in a manner suited for possible verification in condensed matter experiments. To 
further reveal the effect of 9 on the phase transition, we compute in section IV all the critical 
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exponents and verify that the phase transitions remain in the universality class of mean-field 
type to leading order in 1/N. The anyonic character of the quasi-particle is demonstrated in 
section V via an explicit calculation of the relevant propagators and commutators. Finally, 
we summarize and point out possible generalizations in section VI. 



II. RESULTS IN 1+1 DIMENSIONS 



To leading order in 1/N the effective potential for the cr-field is given by the standard 
one-loop formula, in dimensional regularization, 



r ,,2e 2 i oo „ .Jl-2e, , 



(2tt)' 



(5) 



where k n = k n + 2n9 / L and k n = ir(2n + 1)/L. Since an additive constant in 14ff[c, 9] is 
irrelevant, we can first calculate 



dV cS [a,9] 
da 



= Na 



1 



r /i 2e 1 00 



(6) 



After some straightforward manipulation we obtain explicitly 

- ? (- + 1e -1h4tt + ln(/i 2 L 2 )) - C (9) -Y.C 2 M{°LY V , (7) 



dV eS [a,6\ Na\n/j, 2€ 1,1 



da 
where 



7T 



gl 



1 



Co(o) = - £ C(i + y) (1/2 + ey + (1/2 - ey 
z i/=i 



9 ^ ln(n - 



1/2-9 



+ 



1/2 + 9 



2 Ln(n - 1/2 + 0) n(n - 1/2 - 9) 



and 



C 2 , (9) = Li£f M (j^ [^ 2z/ + V2 + *) + C(2^ + 1, 1/2 - 



(8) 



(9) 



In the above equation ^ E is the Euler-Mascheroni constant and ((a, z) = J2^=o( n +z)~ a is the 
generalized Riemann Zeta function. By further choosing the modified minimal subtraction 
scheme, the final result for the renormalized effective potential is then given by 



dV eS [a,9] Na 



da 



7T 



-ln(ML) - lE + hi47r - C (9) - £ C 2u (9)(aL) 



2;/ 



(10) 



where M is a shorthand notation for the standard Aj^g. Since the convergence radius of 
the series in oL in the above equation is 2n, it is sometimes necessary to use the integral 
representation of the series 
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C {9) + Y j C 2u {9)z 2v = ln47r - 7£ -ln^+ y(-l) n cos(27m^) / — exp(-z 2 t 
„=i „.=i J ° t 



-) 



(11) 



in Eq.(0). 

The mass gap is solved from dV e ff(cr = m)/da = 0, yielding 

m(L 9) = I Mex p[^ =1 (-l) n cos(27r^) / °° f exp(- 2 2 t - n 2 /4t)] , when L > L c (0), 
1 ' ' \ 0, when L < L c (0), 



where z = m(L, 6 I )L and L C (6 I ) is the critical length of the ring defined by 



7T 



L c{9) = ^exp 



-7 B + ln4-C (^) 



'121 



(13) 



As defined in Eq.(|8|) Cq(9) is a monotonically increasing function of 9 6 [0, 1/2) starting at 
Cq(0) = In 4 and diverging at 9 = 1/2. Thus the longest critical length results when 9 = 0, 
L c (0) = 7rexp(-7£)/M « 1.7639/M, and the shortest critical length is zero, L c (l/2) = 0. 
For other values of 9 the critical length interpolates continuously between these two extrema. 
In Fig.|l] we plot m(L, 9) as a function of 1/L for several values of 9. 

The non-trivial dependence of the mass gap on 9 and L has an interesting implication 
with respect to the small- volume expansion, proposed by Liischer J7|||, for asymptotically 
free theories with phase transitions. Liischer's approach is based on the observation that 
when the volume is small, the running coupling becomes small due to asymptotic freedom, 
and perturbative calculations can then be applied. The physical result is finally obtained by 
doing a sequence of calculations with increasing volumes and extrapolating to the infinite 
volume limit. However, there are some obstacles in doing the infinite volume extrapolation. 
One is due to the existence of instantons, mostly addressed so far in the literature p. The 
other is related to the existence of phase transitions, such as the deconfmement transition 
in QCD at finite temperature. When a phase transition is present the infinite volume 
extrapolation can hardly be smooth. The 1+1 dimensional Gross-Neveu model in the large- 
iV limit provides a perfect example for illustrating the problem. 

Since M = is the only dimensionful parameter, the transition point L c is of order one 
in units of 1/M when 9 = 0. This transition is of course the same as the finite temperature 
phase transition. From Fig.|l| we see that the entire regime of finite mass gap is located in the 
non-perturbative region, where the perturbative running coupling is very large or divergent, 
as exemplified by the first order (dotted line) and second order (dashed line) perturbative 
results, respectively. When the size of the box L is smaller than L c , chiral symmetry forbids 
the mass gap generation. Hence the small-volume expansion in the most naive form does not 
work at all. However, as we have shown, we can delay the phase transition to an arbitrary 
point by varying the relevant boundary condition. This widens the window in which the 
mass gap curve overlapps with the perturbative regime. More concretely, let us examine the 
mass gap behavior near the bosonic limit. To order (9 — 1/2) 2 the mass gap equation is 



ln(^) + ln4.- 7 ^(4(,-l/2) 2 + I^)-- 



(14) 
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which yields, in the region of \9 — 1/2| < g 2 (L)/27r 

9 2 (L) 



m(L,6) 



An 2 (9-l/2) 2 
L \ 7{L) 



(15) 



where g 2 {L) = it/ ln(47re _7£; /ML) is the perturbative running coupling constant. The crit- 
ical region is then shifted to g 2 {L) ~ \6 — 1/2 1. Thus, close to the bosonic limit 9 = 1/2, 
the mass gap curve possesses a big region in which a perturbative calculation can be carried 
out. At the very value of 9 = 1/2 the gap equation Eq . (|I"0|) becomes identical to that of the 
non- linear a- model in the large- N limit [[?]]. In practice, some intermediate values of 9 may 
be optimal to make the curve flat enough for a better extrapolation, though the choice of 9 
would only be known a posteriori. Using more general boundary conditions, such as twisted 
boundary conditions [ID], to improve the convergence in the Yang-Mills gauge theory have 



been pursued JTI] . 



Strictly speaking, the analogy with finite temperature phase transitions is no longer 
valid when more than one spatial directions are compactified. However, even in these cases, 
rapid cross-overs are still expected, though there would be no true singularities, which only 
associate with phase transitions. 

We would like to mention in passing that the phase transition in the 1+1 dimensional 
Gross- Neveu model with a finite L is no longer possible when N is finite due to the existence 
of the a- kinks, as shown in Ref . [[H| . For this reason we will not dwell further into details on 
the thermodynamical observables in the 1+1 dimensional case. However, since this absence 
of the phase transition is invisible in the large- iV expansion to any finite order, the results 
we found above are entirely legitimate within the domain of the 1/N expansion. As we will 
show in the next section, qualitatively similar results are obtained in the 2+1 dimensional 
Gross-Neveu model, where a phase transition can exist for an infinitely long cylinder (but 
with finite circumference L) at zero temperature. 



III. RESULTS IN 2+1 DIMENSIONS 

The effective potential for the a- field in 2+1 dimensions to leading order in 1/N has a 
similar form as Eq.(0) 



V eS [a,9}=N 



2t„2 i oo 



T E 



2c 



d 2 ~ 2e UJ 



(2+ 



2-2e 



\n(u 2 + k 2 n + a 2 ) 



(16) 



Again we first calculate the derivative of V e s with respect to a. After some standard ma- 
nipulation we obtain 



where 



da 



C 2v {9) 



Na 



r U 2e 1 

t-Y + ^—f ln(4 cos 2 7r6 



9% 2ttL 



2-kL 



2;/ 



-1 



v ^2tt 



1 \ 2v 



C(2v,l/2 + 0) + C(2v,l/2 



(17) 



(18) 
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Due to the use of dimensional regularization we do not encounter infinities here. However, in 
order for the theory to possess spontaneous parity breaking in the limit L — > oo we demand, 
as in Ref. pL 



2c 



9% 



M 
2^ 



(19) 



M will turn out to be the dynamical mass of the particle in the limit L —>■ oo. We want 
to remind the reader that in 2+1 dimensions g 2 B has the dimension of inverse mass. If we 
further use the infinite product representation of 



cosh z — cos 2na = 2 sin 2 ira (l + 



An 2 a 2 



no 



fc=i 



4ir 2 (k + a) 



1 



An 2 (k - a) 2 



(20) 



the series in Eq.flTTD can be summed into a closed form 



dV eS [a,9] Na 



M + - ln(2 cosh aL + 2 cos 2n6 
L 



da 2tt 

The mass gap is immediately solved from 8V c ^(a = m)/da = 



(21) 



m{L,6) 



j cosh 1 
0, 



1 + \{e ML - e ML ^))\ , when L > L c {6), 

when L < L c {6), 



(22) 



where the critical length L c {9) is defined by 



L c (9) 



0, 

j-; ln(4 cos 2 7r6 



when 9 e (1/3,2/3), 
otherwise. 



(23) 



The free energy, equal to the effective potential evaluated at a = m(L, 9), has the form 



F[L,9] 



N 
2^ 



M 



m 2 (L,9) + 



m(L,8) 



da a ln(2 cosh aL + 2 cos 2n9) 



(24) 



For convenience from the experimental point of view we would like to regard 9, rather 
than l/L, as the independent variable, because one could use the same sample and only 
change the magnetic field strength in the solenoid. In Fig.0 we display the mass gap as a 
function of 9 for three values of l/L. The distinct feature in the figure is the non-analyticity of 
m(L, 9) as a function of 9 when L < L c (0), reflecting the fact that the Aharonov-Bohm effect 
pushes the critical length of the phase transition to a smaller value or effectively changes 
the mass scale of the theory to a larger value. In Fig|| we plot the phase diagram in the 
(l/L, #)-plane determined from Eq. (|22|) . We not only find the change of the critical length 
as 9 is varied, we also find, surprisingly, that there is a region 9 6 (1/3,2/3) where parity 
is always broken, no matter how small L is. In other words, even when the renormalized 
mass scale M = 0, the system would still prefer to stay in the parity breaking phase for 
9 e (1/3,2/3). 

To understand why there should be such a region of 9 where parity is always broken, let 



us focus on the cx-propagator. We know from the literature El that the mechanism of parity 
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breaking at zero temperature (corresponding to L — > oo and 9 = in our case) is due to the 
formation of a quasi-bound state, which shows up as a "pole" at the fermion-antifermion 
threshold, u 2 = AM 2 . It is this quasi-bound state, analogous to the Cooper-pair in an usual 
superconductor, that condenses and renders the expectation value (^0) non-vanishing. By 
an explicit calculation of the a-propagator, we can show that this quasi-bound state persists 
whenever the mass gap m(L, 9) is finite. Therefore the mechanism for parity breaking is 
the same when 9 e (1/3,2/3) as when L —>■ oo and 9 = 0. Since the threshold becomes 
zero in the symmetric phase due to the vanishing mass gap, all states have to lie within the 
continuum and hence there can not be any bound state. One can easily verify that uj 2 = 
is no longer a pole for the a-propagator in the symmetric phase. 

Qualitatively, the symmetric phase is a kinetic energy dominated phase, while the sym- 
metry breaking phase is the interaction dominated phase. When the size of the cylinder L 
becomes smaller, the kinetic energy becomes more and more important, if none of the dis- 
crete momentum levels vanish. In our case the effective momentum is k n = n(2n + l + 6)/L, 
which interpolates between the fermionic and bosonic cases. When 9 is finite, one of the k n 
is closer to zero and thus the role of the kinetic energy is reduced. This in turn implies that, 
effectively, the interaction becomes stronger when 9 increases from the fermionic value to 
the bosonic value 1/2. We will show, in section V, that the presence of 9 indeed modifies 
the phase acquired when interchanging two particles, from the original fermionic value, it, 
to 7r + 2tt9. This extra phase could be interpreted as an additional contact interaction. One 
such example in the context of non-relativistic quantum mechanics was noticed some time 



ago by Leinaas and Myrheim |14 



Another physical observable is the induced (or persistent) surface current density circling 
around the cylinder, defined by 

_ dF[L, 9] L dF[L,6] _ N siu2tt9 r™(L,8)L x 
^ ' ^ = dai ~ ~2~ir 39 2ttL 2 Jo % cosh x + cos 2tt9 ' ^ ' 

One can easily recognize that the overall factor sin 2tt9 in the above equation has the same 
origin as a similar factor in the Aharonov-Bohm scattering amplitude || . This overall factor 
implies that the induced current vanishes at the fermion and boson limits or when 9 is integer 
or half integer, in analogy to the vanishing of the Aharonov-Bohm scattering cross section in 
the same limits. A similar induced current on a ring was derived in terms of a single particle 
picture in Ref. |l5j more than a decade ago. What is novel in our case is that m(L, 9) 



is a dynamically determined quantity, due to many-body effects, rather than an external 
parameter. In Fig.^ we plot J[L,9] as a function of 9 for several values of l/L. Again we 
observe the non-analyticity, inherent from the non-analyticity of m(L,9). 

Expanding the free energy in Eq.(E¥) as a power series in m(L, 9)L, we obtain 



= Nm(L,9) 2 \_ M + 1 cog2 + 1 ~ c / , 2u ^ ( 2 . ) 

An L L , ( v ' 



with 



■ i - 2u 



CV(0) = V^nOrJ W p ' V2 + 0) + C(2u, 1/2 - 9) . (27) 
v[y + 1) x 2tt / 1 J 
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This expression for F[L, 9} can be interpreted as a virial expansion, in analogy to the expan- 
sion of pressure as a power series in the particle density. The corresponding virial coefficients 
are C2 V {9\ which become singular when 9 = 1/2, since ((a, z) ~ z~ a when z is small enough. 

We will show in section V that 9 plays the role of the statistics parameter, with 9 = 
being the fermionic limit and 9 = 1/2 being the bosonic limit. The divergence of the virial 
coefficients at 9 = 1/2 apparently reflects the fact that it is not possible to use a single series 
to represent the free energy with all the coefficients smoothly interpolating between the 
fermionic and bosonic regimes, even though both regimes are well behaved by themselves. 
Put differently, the fermionic regime and bosonic regime are separated by non-perturbative 
physics (in terms of the statistics parameter). A much milder singularity (cusp) is found 
for the second virial coefficient in the non-relativistic anyon system considered by Arovas, 



Schrieffer, Wilczek and Zee ]TJJ. One may speculate that the stronger singularity we find is 
rooted in the deep connection between statistics and spin in fully relativistic theories. It is 
worth emphasizing again that the non-analyticity is a manifestation of the Aharonov-Bohm 
effect in an infinite many-particle system. Without a phase transition, the 9 dependence of 
physical observables would never be singular. 



IV. UNIVERSALITY 

In the last two sections we found that many thermodynamical observables are strongly 
influenced by the Aharonov-Bohm effect. In order to reveal the precise influence of the 
Aharonov-Bohm effect on the phase transition we calculate the critical exponents in this 
section. It turns out that if we fix 9 and vary 1/L, in analogy to the usual variation 
of temperature, the universalities of the chiral restoration in the 1+1 dimensional case 
and parity restoration in the 2+1 dimensional case are not altered by the presence of 9 
to leading order in l/N, as long as the phase transition is not totally destroyed. Since 
the calculation process is insensitive to the space-time dimensions involved, we will only 
present the calculation in the case of 2+1 dimensions. For notational simplicity let us define 
r = [L — L c (9)]/L c (9), with L c {9) given by Eq . (|23|) . For detailed definitions of the critical 
exponents we refer to the book |I7| |. Because the phase transition exists when 9 E [0, 1/3) 



and 9 E (2/3, 1], the critical exponents are defined only in these regions in 2+1 dimensions. 
From Eq.(|22"D we can immediately solve for the mass gap m(L, 9), for small r > 0, 

m{L, 9) « M r 1 / 2 , (28) 

'ln(4 cos 2 ti9) 



which in turn yields (3 = 1/2. 

The free energy clearly vanishes in the symmetric phase when L < L c (9), since m(L, 9) = 
there. In the symmetry breaking phase the free energy is given, as can be seen directly 
from Eq.pj), by 



4tt 



M + ^p-M + C 2 {9) (m{L, 9)L 



T 



2 



NM 3 1-C 2 (0) 
47r ln(4 cos 2 n9) 

(29) 
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with 62(0) = 1/8 cos 2 7T0. The specific heat c v is proportional to the second derivative of 
F[L, 9} with respect to r. Thus we see that c v is finite in both phases and has a discontinuity 
at the critical point, which implies that a = a' = 0. 

To obtain the exponents 7, v and 77 we need to calculate the two-point correlation func- 
tion for the ex-field in momentum space in the "static" limit, D a (u = 0,p). Here u should 
be understood as the momentum in the x-direction and p represents the momenta along the 
Euclidean time and y-direction, respectively. To leading order in 1/JV this correlation func- 
tion is given by simple iteration of the bubble graph. After some lengthy but straightforward 
calculation we obtain, to order p 2 , 

D.( U = 0, p) = £ [ g^W pS + aM r] . (30) 

where a = 1 when r > and a = —1/2 when r < 0. At the criticality, or r = 0, 
D a {ui = 0,p) oc p~ 2 , which implies r\ = 0. The susceptibility x is proportional to D a {uj = 
0,p = 0) oc t~ 1 in both phases. Thus, we have 7 = 7' = 1. Eq.(|30D also implies that the 
correlation length £ has the form 



or equivalently v = v' = 1/2. 

In order to obtain the last exponent 5 we need to introduce an uniform external field h 
coupled linearly to ipip (or a mass term). In the presence of h the gap equation to leading 
order in 1/JV is still given by Eq.flTTD except for the replacement a — > a + h in the one- loop 
part. As a consequence of this replacement the gap equation now becomes 

{1 00 21 
^mW)^ 02 ^^ 1,9 ^ ) h ' (32) 

At the critical point we have m 3 oc h, or in other words, 6 = 3. 

This completes the calculation of the exponents in 2+1 dimensions. The same result was 
obtained in Ref. |18| for 6 = 0. One obtains identical exponents by a similar calculation in 
1+1 dimensions. We summarize the final result in Table |I[ It hardly escapes notice that 
the exponents in Table [I] are of the mean-field type. This does not constitute a surprise 
since we only work up to the leading order in 1/JV. However, one should be aware that 
the large- JV limit alone does not guarantee a mean-field result. For example, the non-linear 
cr-model in lower than 4 dimensions violates mean-field universality even in leading order in 
1/JV ||17|| . The real reason in our case is that no <r-loop contribution is involved to leading 
order in 1/JV. Of course, we anticipate that the degeneracy of the critical exponents in 1+1 
and 2+1 dimensions is lifted when we include higher order corrections. It perhaps should be 
mentioned explicitly that, though the critical exponents are independent of 9, the amplitudes 
do depend on 9. 

Finally, we would like to point out that, if we regard 9 as the independent variable and 
fix 1/L, we find that the functional dependence on 9 of some physical observables, such as 
the mass gap and the induced current, near the critical region depends on the value of 1/L, 
as explicitly indicated in Fig.0 and Fig.HI, with L = L c (0) as a special point. 



10 



V. THE NATURE OF THE QUASI-PARTICLE 

Since the interactions between quasi-particles are of the order 1/N, which is vanishingly 
small in the large N limit, we are allowed to discard the interaction all together when we 
attempt to elucidate the nature of the quasi-particle. The effective Lagrangian is simply 
given by 

£ eff = i)[id^ - a^flJfr'V - rm/jip , (33) 

with m = m(L,9) and ij)(x + L) = —ip(x). The constant gauge field can be eliminated 
by a gauge transformation, as mentioned in the Introduction, at the expense of the ^-field 
becoming multi- valued when 9 ^ or 1/2. This multi-valuedness signals that the quasi- 
particle is anyon-like. 

To be more specific, let us calculate the propagator for the ^-field, which can be easily 
written down in momentum space as 

S{k n ,uj;6) =%— r -. (34) 

uj z — K„ — m z + it 

The propagator in coordinate space is then obtained by a Fourier transform 

S{x, t-6) = - J2 T~ e t(uJt - knX) S(k n , cu; 9) = [{id, - a,{9))^ + m)S (x, t; 9) , (35) 

L „_ ' J-oo Z7T 



n=— oo 



where Sq(x, t;9) is the "bosonic" propagator 

S Q (x,t;9) = - V e~ iknX f°° —e wt ^ , (36) 

L n= _oo •'-oo 2% uj 2 — k 2 — m 2 + it 



n=—oo 
oo 



1 °° 1 r i 

= T E e- iknX — \9{t)e- lEnt + 9(-t)e iEnt ] , (37) 

L n=-oo 2E n 1 J 

and E n = \fk 2 + m 2 . To carry out the discrete momentum sum we use the Poisson summa- 
tion formula 

OO OO /*OG 

E /(*) = E / tye-nt+fffl. (38) 

n=-oo l=-oo 00 



It is then straightforward, following the procedure in [19|], to obtain 



S (x, t-9)= ]T e inl e i27r(l+x/L)e S B (x + IL, t) , (39) 



l=— oo 



where S B (x,t) is given by 

S B {x,t) = ^9{x 2 - t 2 )K (mVx 2 -t 2 ) - i-9{t 2 - x 2 )H {2 \m\/t 2 - x 2 ) . (40) 

2 77 4 
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In the above equation t should be understood as t(l — ie) = t — iesign(t) or equivalently 
t 2 — > t 2 — ie, to avoid divergences on the light-cone. 

The physical meaning of each term in Sq(x, t; 9) can be easily recognized. The "/" sum is 
obviously a sum over winding numbers. The first factor is due to the anti-periodic boundary 
condition. The second factor accounts for the total flux embraced by a quasi-particle with a 
path which winds around the origin times. The third factor is nothing but the relativistic 
boson propagator in flat 1+1 dimensional Minkowski space from point Xi = (0,0) to point 
xj = (x + lL,t). Therefore, the above equation is a decomposition of S (x,t; 9) into a sum 
over paths of distinct homotopy class in the universal covering space with 9 taking the role 
of the statistics parameter. In the limit L » t all terms except I = are exponentially (or 
algebraically when m — 0) suppressed and S(x, t; 6) approaches the usual boson propagator. 

The interpretation of 9 as the statistics parameter can be also understood from the 
point of view of interchanging a pair of quasi-particles. Since these quasi-particles can not 
penetrate through each other due to the hard core (inherent from the fermionic Lagrangian 
when 9 = 0), the only way to interchange the positions of quasi-particle number 1 located 
at x\ = and quasi-particle number 2 at X2 = x is to let quasi-particle 1 travel through the 
interval (0, x), while quasi-particle 2 travels through the interval (x, L) and then (L := 0, x). 
The combined world line exactly circumscribes the ring once and therefore accumulates a 
phase factor 

expji j dx^(9)\ = e i2wd . (41) 

In possession of the explicit expression for S(x, t; 9) we can use the Bjorken- Johnson-Low 
formula f2(| to calculate the equal-time anti-commutation relation 



{^(x,0),^(0,0)}= (lim i S{x,t;9) l0 - lim S(x,t;9) lo ) (42) 

and see whether the usual anti-commutator is affected by the presence of the magnetic flux. 
It is readily verified that 



oo 

limlimS^a;, t(l — ie); 9) = -sign(t) V] ( — + IL)^q + terms independent of t , (43) 

e— >0 1— >0 2 , — ' 

t=— oo 



which in turn yields 

oo 

{Mx, 0), 4(0, 0)} = 5 aP £ {-l) l 6(x + IL) . (44) 

£=— oo 

While the sum over "/" is again due to the geometry of the ring, we have verified that 9 does 
not enter the canonical equal-time anti-commutator in our model, in contrast to the graded 
commutator obeyed by anyons induced via the Chern-Simons action. This result should not 
be surprising in light of the local character of the canonical commutation relation and the 
fact that the magnetic flux is global and never attached to each particle in our model. On 
the other hand, the propagator contains global information and should have a non-trivial 
^-dependence, which in turn yields the non-trivial dependence on 9 of thermodynamical 
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quantities. Thus, the nature of the quasi-particle now becomes clear. Locally, the quasi- 
particle is a fermion, but globally, the quasi-particle behaves like an anyon. 

Similarly, the propagator in 2+1 dimensions can be calculated and is again given by 
Eq . (|3~5l) except that S (x, t; 9) in Eq . (|39|) should be replaced by the 2+1 dimensional version 

oo 

S (x, y, t;9) = £ e M e^ l+x / L ^S B (x + IL, y, t) , (45) 



l=— oo 



with 



f> —m\/r' 2 —t' 2 p—imsfW—r 2, 
■"'4.2 2\ e 



S B (x,y,t;6) = 9(r 2 -t 2 ) - i9(t 2 - r 2 ) (46) 

and r = ^x 2 + y 2 . Using the identity 

i 1 ^ (r2 + 6 e 2 )3 /2 = ^(x)6(y) , (47) 
the anti-commutator in 2+1 dimensions is then given by 

oo 

{^(x, y, t), 4(0, 0,t)} = 5 aP (-!)^0 + lL)6(y) . (48) 

l=— oo 

It is not difficult to see that the quasi-particle in the 2+1 dimensional case is anyon-like only 
in the x direction. In other words, only those world lines which circumscribe the axis of the 
cylinder non-trivially can acquire phases of integer multiples of 2tt6. 

VI. SUMMARY AND OUTLOOK 

We considered the statistical mechanics of the Gross-Neveu model on a ring and on a 
cylinder with a magnetic solenoid coinciding with the axis. Among the interesting results 
we obtained are 1) the periodic but non-analytic dependence of thermodynamical observ- 
ables on the magnetic flux (9) and 2) the existence of an interval of 9 G (1/3, 2/3) (modulo 
integers) where parity is always spontaneously broken. All these phenomena are explicit 
manifestations of the Aharonov-Bohm effect in thermodynamics or in many-body systems. 
We further showed that the mean-field nature of the existing phase transitions is preserved 
to leading order in 1/N, by verifying the ^-independence of all the critical exponents. The 
precise nature of the quasi-particle, locally fermion-like and globally anyon-like, was illumi- 
nated through the calculation of the equal-time commutator and the decomposition of the 
propagator into a sum over paths of distinct homotopy class or winding number. 

There are several directions in which our present work can be extended within the same 
model. The most interesting one is to include an external space-time dependent electromag- 
netic field and study the response to it. To set up junctions and investigate tunnelings and 
interferences could also be interesting. Although we only worked to leading order in 1/N, we 
do expect that our results, especially in the 2+1 dimensional case, are qualitatively stable 
against higher order 1/A-corrections, based on the next to leading order study of the same 



13 



model in Ref. [£1]] in the absence of 9. Nevertheless, an explicit calculation of the higher 
order corrections should be carried out. One may also try to introduce a chemical potential 
and study the situation away from the half-filling point. 

As we have been emphasizing all along, the periodic but non-analytic ^-dependence 
of thermodynamical observables is a manifestation of the Aharonov-Bohm effect in many- 
body systems with phase transitions, not necessarily pertinent only to the Gross-Neveu 
model itself. For this reason we anticipate that similar phenomena can be found in other 
models. The detailed microscopic interactions are likely to play minor role, as long as phase 
transitions exist. Thus, in order to make direct contact with experimental verification, it 
would be very interesting to consider more realistic models. 



14 



REFERENCES 



[1] R. E. Prange and S. M. Girvin, The Quantum Hall Effect (Springer, New York, 1987); 
A. Shapere and F. Wilczek, Geometric Phases in Physics (World Scientific, Singapore, 
1989) 

[2] F. Wilczek, Fractional Statistics and Anyon Superconductivity (World Scientific, Sin- 
gapore, 1990) 

[3] Y. Aharonov and D. Bohm, Phys. Rev. 115 (1959) 485 
[4] D. J. Gross and A. Neveu, Phys. Rev. D10 (1974) 3235 

[5] E. Fradkin, Field Theories of Condensed Matter Systems (Addison- Wesley, Redwood 

City, 1991) chapter 4 
[6] X. G. Wen, F. Wilczek and A. Zee, Phys. Rev. B39 (1989) 11413 
[7] M. Liischer, Phys. Letters B118 (1982) 391 
[8] M. Liischer, Nucl. Phys. B219 (1983) 233 

[9] P. van Baal and A. S. Kronfeld, Nucl. Phys. B (Proc. Suppl.) 9 (1989) 227; 
H. Hansson, P. van Baal and I. Zahed, Nucl. Phys. B289 (1987) 628; 
J. Roller and P. van Baal, Nucl. Phys. B302 (1988) 1 
[10] G. 't Hooft, Nucl. Phys. B153 (1979) 141 

[11] A. Coste, A. Gonzalez- Arroyo, C. P. Korthals-Altes, B. Soderberg and A. Tarancon, 
Nucl. Phys. B287 (1987) 569; 

M. Garcia Perez, A. Gonzalez- Arroyo and P. Martinez, preprints INLO-PUB-18/93, 
FTUAM-93/40 (1993) 
[12] R. F. Dashen, S. K. Ma and R. Rajaraman, Phys. Rev. Dll (1975) 1499 
[13] B. Rosenstein, B. J. Warr and S. H. Park, Phys. Rep. 205 (1991) 59 
[14] J. M. Leinaas and J. Myrheim, II Nuovo Cimento 37B (1977) 1 
[15] C. C. Gerry and V. A. Singh, II Nuovo Cimento 73B (1983) 161 
[16] D. Arovas, J. R. Schrieffer, F. Wilczek and A. Zee, Nucl. Phys. B251 (1985) 117 
[17] C. Itzykson and J.-M. Drouffe, Statistical Field Theory (Cambridge Univ. Press, Cam- 
bridge, 1989) 

[18] B. Rosenstein, B. J. Warr and S. H. Park, Phys. Rev. D39 (1989) 3088 

[19] N. N. Bogoliubov and D. V. Shirkov, An Introduction to Theory of Quantized Fields 

(Wiley-Interscience, New York, 1959) 
[20] J. D. Bjorken, Phys. Rev. 148 (1966) 1467; 

K. Johnson and F. E. Low, Prog. Theor. Phys. (Kyoto) Supp. 37-38 (1966) 74 
[21] G. Gat, A. Kovner, B.Rosenstein and B. J. Warr, Phys. Lett. B240 (1990) 158 



15 



FIGURES 



FIG. 1. Mass gap m(L, 9) (in units of M) as a function of 1/L for 9 = (fermion limit), 0.3, 0.4 
and 1/2 (boson limit) in the 1+1 dimensional Gross-Neveu model. Also plotted are the first order 
(dotted line) and second order (dashed line) perturbative results. The sudden turns of the dotted 
and dashed lines for small values of 1/L reflect the fact that the perturbative running coupling 
diverges at L = 4L C (0). 

FIG. 2. Mass gap m(L,9) as a function of 6 for L = 1, L = L c (0) = In 4 and L = 2 (in units 
of M). The functional form of the ^-dependence of m(L,9) in relevant regions is also indicated 
explicitly. 

FIG. 3. Phase diagram in (1/L, 0)-plane (with M = 1). In the interval of 9 £ (1/3,2/3) the 
system is always in the parity broken phase, independent of the value of 1/L. 

FIG. 4. Induced current J[L, 9] as a function of 9 for L = 1, L = L c (0) = In 4 and L = 2 (with 
M = 1). The functional form of the ^-dependence of J[L,9] in relevant regions is also indicated 
explicitly. 
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TABLES 



TABLE I. critical exponents 
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